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I. INTRODUCTION 

The interplay of disorder and quantum criticality is a 
long-standing and still open problem in condensed mat- 
ter theory. One aspect of this problem which has re- 
ceived considerable recent attention is the 'quantum Grif- 
fiths' behavior which has been shown to occur |H|ea|t,quan- 
tum critical points in certain model systemsH'clclaa. The 
model systems in which quantum Griffiths behavior has 
been unambiguously demonstrated all possess a crucial 
common feature, namely that in the absence of disor- 
der the critical degrees of freedom exhibit dissipationlcss, 
Hamiltonian spin dynamics (indeed typically character- 
ized by dynamical exponent z — 1). However, many sys- 
tems of experimental importance involve magnetic, jde.- 
grees of freedom coupled to conduction electronaju^BEl, 
and therefore overdamped dynamics implying a pure- 
system critical behavior characterized by z > 1. Ex- 
tension of the theory of quantum Griffiths behavior to 
this c ase , is . t herefore an important issue. In a series of 
paperst2lOll3 Castro-Neto and Jones have argued from 
various points of view that such overdamped systems ex- 
hibit quantum Griffiths behavior similar to that exhibited 
by undamped systems, and they and others have further 
argued that this phenomenon is at the heart of the 'non- 
fermi-liqiiidl, behavior observed in many heavy fermion 
materialsoO. 

In this paper we examine the issue of quantum Grif- 
fiths behavior in nearly critical systems exhibiting over- 
damped dynamics, finding that it is essentially nonex- 
istent, being replaced instead by 'superparamagnetic' 
behavior. The essence of our analysis is this: in un- 
damped models quantum Griffiths effects arise from an 
interplay between the low probability of nucleating mag- 
netic 'droplets' in the paramagnetic state and a low but 
non-negligible quantum tunnelling of these droplets. In 
a metallic, dissipative environment there is a strong sup- 
pression of tunnelling by dissipation, so that the droplets 
which dominate the susceptibility behave more or less 



classically, leading to superparamagnetic behavior rather 
than quantum Griffiths behavior. 

Our cesults amount to an implementation of ideas out- 
lined inEj and to a generalization, to a non-vanisliing den- 
sity of defects, of a previously reported analysiaij of the 
'magnetic droplet' produced by a single, spatially local- 
ized defect, and rely heavily on the results of this previous 
work. The method used to analyse the dynamics of a dis- 
tribution of defects is similar to that used intil and the 
broad qualitative features of the results we obtain are 
very similar to those obtained in that work. However, 
the specifics and the physical implications seem different. 
The issue is discussed in more detail in the conclusion. 

The outline of this paper is as follows. In section II we 
present the model and the method used in our analysis. 
In section III we show that the approach reproduces re- 
sults previously obtained in the dissipationlcss case. Sec- 
tion IV presents our new results concerning Griffiths-like 
behavior in systems with overdamped dynamics. Section 
V is a summary, comparison to other work, and conclu- 
sion, and is written so that readers uninterested in the 
details of the derivations may obtain from it the essence 
of our results. 



II. MODEL AND METHOD OF SOLUTION 
A. Model 

The canonical quantum Griffiths problem concerns the 
effect of weak disorder added to a 'pure' (non-disordered) 
system which possesses an Ising symmetry and is tuned 
to be near a quantum critical point. We consider a sys- 
tem in imaginary time and 3 spatial dimensions (differ- 
ences occurring for two spatial dimensions warrant a sep- 
arate treatment, which will be presented elsewhere). The 
model is described by the action 

S — S'static + 'S'dyn + "Sdisordor (1) 



2 



with 



^static — — / At 



so 



+fo'[V0(x,T)]V-(/.^(a;,r) 



(2) 



Here </) is a dimensionless scalar order parameter, is 
the basic energy scale of the theory (perhaps of the order 
of the mean Kondo temperature for a heavy fermion sys- 
tem) and is fixed by normalizing the coefficient of the 0'* 
term to unity, ^0 is the basic length scale (typically of the 
order of a lattice constant) , ^ is the magnetic correlation 
length, and /3 is the inverse temperature. It is convenient 
to define a parameter r = (^o/C)^ > which measures 
distance from criticality. We consider only parameters 
such that the pure system is in the paramagnetic phase. 

We take the disorder to couple to the square of the 
order parameter via 



dT 



|fF(x)02(x,r) (3) 
so 



and assume it to be Gaussian distributed with correlator 
((...) represents average over configurations of the disor- 
der) 



y 



Co 



(4) 



where the kernel K{u) decays on the scale m ~ 1 and 
satisfies / (fiuK(u) = 1. Because we are interested only 
in length scales x — y > ^0 we will take K to he a, S 
function. The dimensionless quantity Vq parameterizes 
the strength of the disorder. Weak disorder corresponds 
to Vo < 1. 

The dynamic term S'dyn is crucial to the quantum crit- 
icality described by Eq. |l| and to our subsequent discus- 
sions. We consider two cases: (i) dissipationless, z = 1 
dynamics^ as is the usually assumed in studies of quan- 
tum Griffiths behavior, with 



r.(2=l) 
dyn 



Eo 
Stt 



dT 



d^x 
so 



eo 



dT 



(5) 



Here c is a characteristic velocity of the undamped ex- 
citations, such that c/^o is an energy presumably of the 
order of Eq. 

(ii) Hertz antiferromagnet, z = 2 dynamics, correspond- 
ing to the generic antiferromagnetic transition in a fermi 
liquid: 



dyn 



j(z=l) 
dyn 



87r£;o ^ r 



So 



where 



(r, w„) ^ Eo dT( 

JQ 



(r,T) e^' 



(6) 



(7) 



In these conventions the dynamics are dissipative (i.e. 
dominated by the T term) if w < w* = c^/(e^r), and non- 
dissipative at higher frequencies. One expects in most 
systems (and finds for example in a weakly coupled fermi 
liquid or in the slave boson theory of the Kondo lattice) 
that all scales are of roughly the same order, i.e. that 

Eo ~ c/Co ^ r. 



B. Method 

1. Overview 

The dissipative term in Eqj6| corresponds to a long 
ranged interaction in time and renders available numer- 
ical methods prohibitively difficult to apply. To analyze 
the model defined by Eq. |^ we use-simple analytical ar- 
guments modelled on those of RefE. We note that the 
effective dimensionality of the model defined by Eq. |l| is 
dcff — d + z. In this paper we consider only the spatial 
dimension d = 3 so we are concerned only with models 
at and above the upper critical dimension dc = 4, so that 
quantal and thermal fluctuations of the order parameter 
in a fixed disorder configuration can be treated by an 
essentially mean field approximation. The usual fluctu- 
ation analysis which justifies the mean-field approxima- 
tion for doff > dc involves a translation-invariant model 
and fluctuations for which momentum is a good quantum 
number. Here we must deal with fluctuations in a sys- 
tem whose translation, invaciance is broken. These were 
investigated in Refs.tj andlij and were found not to af- 
fect the structure of the static mean field solution when 
dcS > dc (except for some insigniflcant changes in some 
constants) . _ 

As noted for example bycJ, in the presence of the ran- 
dom potential, the crucial feature of the mean fleld solu- 
tion is the presence of droplets: regions in which the order 
parameter is locally non-vanishing. Quantum Griffiths 
effects then arise from dynamical fluctuations of these 
droplets; to study them one must estimate the droplet 
density and tunnelling rate. We use statistical arguments 
and mean field analysis to estimate the density and an 
adaptation to the present case of the analysis presented 
for a droplet produced by a single point defect in Ref.Ej 
to estimate the tunnelling rate. 



2. Probability for the existence of a droplet 

The assumption that the model is at or above its up- 
per critical dimensioH-|ip«ans that mean-field theory is 
a good starting pointtjila. We therefore consider static 
configurations, 4>{x), which minimize the combination of 
Eqs.pjand (||). These satisfy 



^oV^ 



(x) + r^{x) + (l){xf = ^V{x)<j>{x) (8) 



If V{x) — 0, then because we assume r > the mini- 
mum corresponds to 4'{x) = 0; however regions in which 
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V{x) < can lead to 0(a;) 7^ 0. In the regions where 
V{x) = const. < 0, (/){x) is roughly constant whereas 
in between these regions ^(x) decays exponentially. We 
refer to the regions where </> is not exponentially small 
as 'droplets'. If the droplets are reasonably dilute, one 
may set = in the exponential tail rcgional3 and es- 
timate the density of droplets of a given size and mean 
amplitude. 

To motivate our estimate we first consider solving Eq ^ 
if V{xX-r^ V for |x| < R and Vq = otherwise. A previous 
paperE^ considered a special case of this equation, with 
Vo{x) = VS^'^\x) and the solutions found in that work 
may easily be modified for the present case. In d = 3 
one finds that the solution is, roughly (and neglecting 
unimportant logarithmic factors in the x-dependence) 



</.(x) 



for X < R 
for R < X 



(9) 



In other words, the magnetic order induced by the re- 
gion of attractive V is roughly constant inside the region 
and decays outside it, initially as 1/a; and exponentially 
for distances larger than a correlation length from the 
boundary of the attractive potential region. Inserting the 
above ansatz, Eq. |^, into Eqs. ^ and ^, and minimizing 
the resulting action with respect to c/jq yields 



t7 _ £0 



V = ^a(i?/e) + foKR/O 



(10) 



with a{x) = 1 + 3/x + 3/x'^ and b{x) = 1 + 3/x - 
4e'*^r[0, x] . These particular forms for a, b depend on the 
specific potential configuration studied (here V — const 
for X < R and V = otherwise) and on the varia- 
tional approximation used; but we argue that a generic 
droplet is described by a similar equation with a, b func- 
tions which vary on the scale i?/^ ~ 1 and which tend to 
unity as i?/^ 00. Also in 3 dimensions a{x) as 
X ^ while b{x) tends to a constant for x < The 
precise forms of a, b affect only nonuniversal details such 
as widths of crossover regions. In this paper we shall 
assume 



a{x) = 1 -|- 3a; ^ 
b{x) = 1 



(11a) 
(lib) 



where the 3 arises from the difference in integrating a 
constant or 1/r^ over r^dr. 

One sees from Eq ^ that in order to obtain a solution 
at all the average potential, must be smaller than a 
(negative) i?— dependent critical value, 



(12) 



which tends to /^^ as i? ^ 00 and to a number of order 
1 as i? — > ^0- As is evident from these formulae, the nat- 
ural scale of the droplets is the correlation length ^ which 
diverges as the quantum critical point is approached. 



Eqs. |l^ and |2| thus imply that one obtains a droplet in 
a region of linear dimension R only if the average value V 
of the potential in that region is larger than a value of the 
order of Vc{R/(,) {Vc is not an exact estimate because it 
pertains to the idealized disorder configuration discussed 
above). The standard estimate of the probability of a 
region of linear dimension R with mean potential V is 



P{R-',V) 



3 T7, W^o) 



3/2 



exp(-(-) (^) 1 (13) 



This argument 



and we therefore argue that the density A^(i?'^, 0q) of 
droplets of amplitude (fiQ and core size R, must be pro- 

portional to ^ exp I — - 

does not determine the preexponential factors (which in- 
volve, e.g. the issue of whether the region of size R con- 
sidered in Eq. |l^ is part of a larger region which can 
sustain a droplet and numerical factors arising from the 
difference between idealized disorder configuration and 
typical one, which we have absorbed into Vq and Vc )■ Be- 
cause some of our subsequent considerations will require 
an estimate of the preexponential factors, we present the 
following arguments to fix them. 

We begin by making a rough estimate of the fraction 
of sites contained in droplets (i.e. of the fraction of sites 
having a (/)q > 0), as a function of distance from critical- 
ity. As noted above, in principle within mean field the- 
ory 00 is non-vanishing everywhere, but we neglect the 
regions where it is exponentially small, in other words 
we set 00 = in the 'inter-droplet' regions. To perform 
the estimate we coarse-grain the theory to the scale ^. A 
given correlation volume i^"^ will have a non- vanishing (/)q 
if the potential averaged over the droplet volume, I4,voj 
is larger than Vc(0 ~ a(l)Co/C^- From Eq |l^ we see 
that the probability P0 that a given correlation volume 
will have a non- vanishing (j)Q, i.e. a V < Vc{£,), is (recall 
K<0) 



Pri. 



1 



1 — erf 



3/2 



Vn 



(14) 



where erf is the error function. 

Clearly, a picture of independent droplets must break 
down if exceeds the percolation probability Ppoic at 
which the set of correlation volumes with non-vanishing 
00 percolate. Use of Eq. ^ and the estimate for 
three dimensional cubic lattices Ppcrc ~ 0.2 shows that 
percolation will have occurred by the time ^ exceeds 
2.8a{l)^ /Vq . (^pcrc is an underestimate because 



pcrc 



droplets larger than ^ may occur). These estimates also 
show that the natural scale for ^ is Vq'^ and strongly sug- 
gest that the probability that a given site is in a droplet 
(of any size) is a function only of the combination £,Vq. 

We therefore argue that the prefactors in the droplet 
density must be such that the total probability of find- 
ing a site in a droplet, Ptot = / dR^d(j)lR^N{R^ , (j^l) 
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must be a function only of and must be of the order 
of Pporc when ^ is of the order of ^porc- This imphes 



Cvo^o 



exp 



[4>i+v,{R/oy 



(15) 



where the factor of R~^^^ ensures the correct scahng with 
^ and the numerical factor Cvg ~ 11.25 ensures that when 
? = Cporc, we have P = Ppcrc- We emphasize that these 
formulae are phenomenological and must in particular 
break down when ^ approaches ^perc- 

It is convenient to adopt a dimensionless system of 
units in which 

(16) 
(17) 







h = 














^0 



(18) 



for which 



N{y\f)dy^df 



-9/2 



X exp 



{P + <V))' 



dy^df 



(19) 



The factor of expresses the fact that if the proba- 
bility of a given site being in a droplet is a function only 
of u, then the density of droplets must be smaller by an 
extra factor of the typical droplet volume 



3. Tunnelling of the droplet for undamped, z — 1, dynamics 

We now estimate the rate wtun at which a droplet char- 
acterized by the mean amplitude 00 and length scale R 
tunnels in the case of undamped, z = 1 dynamics by per- 
forming a variational instanton calculation using Eqsj2[ 
H and the solution Eq. ||. In the simplest estimate one 
assumes that the droplet maintains its shape while col- 
lapsing and re-forming. To estimate the action associated 
with this process we write the droplet solution as 



(I){x,t) = (j){x)'q{T) 
Substitution into Eqs. ||, pleads to 

5'inst 5'kin 'S'barric 



Skin involves the-iintegral of (9*0) over the droplet and 
as noted in Ref £3 involves the 1 /r 'tail' of the droplet in 
a crucial manner; in contrast, the cost S'barrier of creating 
the instanton does not. One obtains 

^r^=C.„^/Va'(.)/|(g)' (22) 

5'barricr — C^.rric.r f y b'iy) 

X / Eodr{-2,j{r)'+7j{T)^) 



Here Ckin and Cbarrior are nonuniversal constants. Ckin 
involves the square of the ratio Eq/{c/^q) of the basic en- 
ergy scale to the kinetic (or zone boundary magnon) en- 
ergy while Cbarricr IS just a number. In the approximation 
we have employed Ckin — Eq^^/c^ and Cbanicr — 1- The 
functions a' = J x^dx(l)o{x)^ and 6' ~ J x^dx(j)Q{x)* are 
functions with behavior similar to a, b; in our explicit cal- 
culations we set a' = a/3 = l/3 + y~^ and b' = 6/3 =1/3 
for simplicity; again different choices affect only nonuni- 
versal details. 

The action associated with one instanton may now be 
determined by a standard minimization of Eqs 22 2^ and 
is 



(z=l) _ 



3„,3 



SMf'y 



(24) 



For the present model in the present approximation the 
nonuniversal constant 5*1 = Ckin Cbarrior /3and d{y) — 
3-\/a'(j/)5'(y) = \/ o,{y) where the last equality follows 
from our simplifying assumptions a' — 3a and b' — 3&. 
The value of controls the width of the crossover regime 
before the universal behavior is reached, and is linearly 
proportional to Eq^/c 

The tunnelling rate is then given by 



LOoe 



-Si„,t(ro) 



(25) 



Here, luq is an attempt frequency presumably of order Eq 
whose value is beyond the scope of this theory. 

To conclude this section we briefly estimate the action 
associated with a different tunnelling mechanism, namely 
nucleation of a domain wall. For small droplets ('core 
size' R less than ^) the important process was shown to 
be collapse and reformation of the entire drople1ii3. We 
therefore need consider only the case R^ We observe 
that by expanding about the static uniform solution one 
obtains a domain wall with width W ~ 't'o^- The kinetic 
term associated with the domain wall motion therefore 
has one fewer factor of the small quantity 00 ^ f/S, ^ 
VqJ, leading to a larger action and hence a smaller rate, 
in the weak disorder, near criticality limit. We note in 
passing that for 0o ^ 1 the powers of R will be the same 
as we have considered but the extra factor of (f>^^ will 



(20) work in the other direction, favoring domain wall motion. 



(21) 4- Tunnelling of the droplet for overdamped, z = 2, 

dynamics 

For overdamped dynamics two jiimQEtant differences 
occur. First, as shown in Refs.El[E30 the damping 
changes the action associated with a single instanton, 
strongly suppressing the bare tunnelling rate relative to 
that found for undamped dynamics. Essentially, the tun- 
nelling is limited by the droplet 's ability to move through 
a viscous medium rather than by its ability to climb 
(23) over a barrier. Second, and much more important, the 
overdamped dynamics leads to a long-ranged (in time) 
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instanton-instanton interaction, which reduces the tun- 
nelling rate further and indeed drives it to zero if the 
damping .exceeds a critical value, as noted by previous 
authorsBEa. 

To calculate the effects of damping we insert the 
ansatz, Eq. ^ into Eq ^. The new term arising from 
the overdamped dynamics is 



Q - 2 

'-'diss — ^ 



,dr] di] At- 



(26) 



with Tm a 'microscopic' time of the order oioj* — (? j (^qF. 
The net dissipative-coefficient 7 is given for the Hertz 
antifcrromagnet byuJ 



7 



1^ / ^^o(a;)^=c,/22;3a"(y)e/eo 



(27) 



The approximations employed in the previous section 
imply that the nonuniversal constant c-y = Eq/T and 
a" = a' . In a generic system one expects all scales to 
be of roughly the same order, so that in particular c-y is 
expected to be of order unity. 

The estimate of 7 is subject to the important caveat 
that the electron bath which causes the dissipation can 
penetrate the entire droplet. A reasonable estimate of 
the penetration depth, Lp, may be obtained by dividing 
the electron velocity, wp, by the magnitude of the order 
parameter; in rescaled units ip/^ ~ VF/{Eof). We shall 
see below Eq. ^ that the parameters are such that the 
electrons can penetrate the entire droplet. 

We have not been able to solve analytically for the 
instanton; instead we estimate the action by inserting 
the variational ansatz 



dr To 



(28) 



into Eqs 6U23 obtaining S = Skin + Sdiss + Sbarrier with 



As previously remarked, we expect the ratios of the 
various dimensional parameters to be of the order of 
unity; also as we shall see below, in this problem the 
important droplets have / ~ so that provided T is 

less than a number of the order unity times ~ f '^/^ 
(within our approximations the precise numerical factor 



(33) 



is vl5) the Tq term is negligible and one has 



To 



and thus 



(34) 



where C2 (=2.283... in the present approximations) is a 
numerical factor of the order of unity arising from com- 
bining the factors in Eqs |2^ and pO| ). 

We observe that for the value of tq given in Eq |33l the 
term written in Eq is larger than iSbarricr (Eq p3|) by 
two powers of the correlation length (provided that the 
quantity / is of order unity or less, as is the case for the 
situations considered here.) Thus, in the metallic case 
and near to criticality, the difficulty in tunnelling arises 
from moving through the viscous medium, not clpibing 
over the barrier. This result was noted previouslyll3. 

The bare tunnelling amplitude is thus 



{z=2) _ 



ujoe 



-sy 



(35) 



and is much smaller than in the dissipationless case, be- 
cause of the factor in the argument of the exponential. 

The st anda ed macroscopic quantum tunnelling 
argumentaiillij'tZI imply that the instanton-instanton 
interaction renormalizes the bare tunnelling rate so that 
if 7 < 1 then the T = tunnelling rate is 



Skin — 



2CMniy^a'iy)f 

-Bo To Co 



(29) 



Sdiss = 2cj-^y^a'{y)f'^ln{cdTo/Tni) (30) 



Sba 



^barrier 
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(31) 



where ln{cd) = J_^^^l/2dxdyln{l -|- (a; — y)^) ~ 
0.1152.... 

Minimization over the instanton duration then leads 

to 



C, ^^^^ ^ Carrier b' ^ (^^^ / ^EoTo f (32) 



Ck 



l5Ck^n 



V ^0 



(36) 



whereas if 7 > 1 tunnelling stops at T = 0. We see 
from Eq. ^ that 7 is a strong function of the droplet 
size and amplitude; droplets which may tunnel (i.e have 
7 < 1) have a very weak amplitude even in rescaled units: 

Eq ^ is a zero temperature result. At T > the 
'Caldeira-Leggett' renormalization is temperature depen- 
dent. The key question for this paper is the temperature 
at which ujtun{T) < T. If 7 > 1 then uitim 

(T) < T at 

all T < Eq, implying that the droplet behaves classically 
at all T. If 7 < 1 then the usual arguments shows that 
^tun{T) drops below T when T becomes greater than 
t^tun{T = 0), so that Eq |3^ gives the temperature scale 
separating a high-T region, in which the droplet behaves 
classically, from the low-T region, in which it behaves 
quantum mechanically. 
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III. ESTIMATE OF QUANTUM GRIFFITHS 
BEHAVIOR 

A. Overview 

The standard Griffiths estimate is that a droplet of 
magnetic moment Md — / d'^r e*'^' 0o(r) {Q is the 
ordering vector) and tunnelhng frequency wtun [R, 4>q\ 
gives rise to a susceptibility Xd proportional to 
Mj/ (wtun + T). The susceptibility of a system with a 
distribution of droplets is then given by 



X(T) - j d^Rd^^o 



N{R:\4>l)Mj[R, M 



T 



(37) 



For a droplet in an antiferromagnetic system, we find Md 
is a random function with magnitude (jy^R-tYve term pro- 
portional to R comes from the boundary of the droplet, 
where the order parameter amplitude is dropping and the 
cancellation over one unit cell of the antiferromagnetic or- 
der is not complete. A different dependence would change 
prefactors but not affect our results crucially. 

It is convenient to introduce an explicit integral over 
frequency, writing 



x{T) = r 



dhj 



UJ - 



T 



(38) 



so that after conversion to dimensionless units we have 



I{oj)= dy'df\eN{y',n)fY5i^-^t 



(39) 

The prefactor ^ ^ in x arises because each droplet has 
magnetic moment of the order of unity and the density 
of droplets is The quantity ^^N has no explicit 

dependence on ^ (see Eq. p^ . 

We will use the delta function to eliminate the / in- 
tegral in / and perform the integration over y either nu- 
merically or via an extremal value argument. 



B. z=l 



Using Eq. GR yields 



(/ N \ 1/3 
ln(i^)\ 1 



(40) 



Substituting this result into Eqj39| yields 

21ni/=^(c.oM r eN{y',P{oj,y))dy 



I{u;) = 



CO Si 



4/3 



d(y)4/3 



(41) 



where N(y, /(w, y)) is N{y, /) (Eq [T^) with / given by 
Eq||. 



In the limit of very low frequency one may use asymp- 
totic methods to analyse the integral in Eq the ex- 
tremum is at 



ymax 

Substitution leads to 
X(T) 



with (restoring units) 



^asymp 



(0 
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3V35i ei^o' 



(42) 



(43) 



(44) 



This is the familiar quantum Griffiths result: if one 
is sufficiently close to the pure system critical point 
(c?(^) < 1) then the susceptibility diverges, with degree 
of divergence characterized by an exponent which ap- 
proaches unity proportional to one power of the inverse 
correlation length. 

Note that the prefactor in Eq ^ rapidly vanishes as 
criticality is approached, so although the susceptibility 
diverges more strongly, the amplitude of the divergence 
decreases. Note further that in the asymptotic limit, 
/ w 1 so that the mean order parameter density (in- 
tegrated order parameter divided by droplet volume) is 
of the order of Thus the picture that emerges is of 
large, weak droplets. 




0.1 
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FIG. 1: Solid lines: calculated effective Griffiths exponent 
for undamped {z = 1) case defined in Eq. ^ plotted vs 
correlation length normalized to correlation length ^perc at 
which droplets percolate , with (from top to bottom) non- 
universal coefficient 5*1 = .1, .3, 1 and frequency uj — W~^ujo- 
Dashed lines: asymptotic result (Eqji^) for same parameters. 



We have evaluated 



des{uj) = 1 



din [/(w)] 
c?ln(a;) 



(45) 
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via a numerical computation of Eq.^. Fig. 1 shows 
doff Lj) as a function of f for w = O.OOlwg and several 
different values of the non-universal parameter 5*1 (solid 
lines) along with the asymptotic limit estimates from Eq 
We observe that for these low frequencies and not too 
long ^ the asymptotic limit provides a reasonable (but 
not perfect) estimate of the effective exponent: relative 
corrections are of the order of {£,/ln{u)o/(^)'^^^)- We see 
also that depending on the value of the non-universal pa- 
rameter 5*1 , the effective exponent may remain above the 
critical value of unity (corresponding to a non-divergent 
susceptibility) until ^ becomes of the order of £,perc- For 
^ of the order of ^perc the standard quantum griffiths 
approximation (independent droplets) breaks down, and 
one must deal instead with the critical singularities ap- 
propriate to a phase transition in a disordered system; in 
other words with the still unsolved problem of the mixing 
of quantum critical and quantum griffiths singularities. 



z=2 



For overdamped dynamics, some droplets will have 7 > 
1 and therefore will not tunnel at all at T = 0. The 
function I{uj) will thus have a contribution proportional 
to 6{uj) leading to the 1/T behavior expected of classical 
droplets. For those droplets which do tunnel we must use 
Eq |3| in Eq M We write 



I{uj) = IoS{u;) + /rest(w) 
with Irest given by Eq B9^ and Iq by 



(46) 



lo = / dy^dffy^N{y\ f )Q{^{y, /) - 1) (47) 



From Eq 38 we see that if /q is appreciable, then x ~ 1 /T: 



this is the superparamagnetism expected from essentially 
classical droplets. 

We begin by estimating /q. The Q function limits the 
/ integration to 



/' > /mi„(2/) 



Co 



y^y^a{y) 



(48) 



Note that for large /min o,{y)- Further, the typical 
scale for / is so that the penetration depth Lp 

of electrons into the droplet is large: ip/C ~ S,^^^ so the 
assumption that electrons penetrate the droplet is indeed 
self consistent. 

Use of Eq H in Eq P gives 



m) = 



30F 




y '^dy 



(49) 




FIG. 2: Ratio of density of magnetization of non-tunnelling 
droplets Jo (Eq|47|) to total density of droplets Itot ~ J duoI{u)) 
for overdamped case and non- universal constant — 0.1 
(larger values of lead to an lo/Itot ~ 1 even for much 
smaller values of ^, as a function of correlation length (not 
normalized to disorder strength) for dimensionless disorder 
strength Vo = 1, (top curve) .7, .5, .3. Note that for all reason- 
able parameters a non-negligible fraction of droplets do not 
tunnel at all. 



Iq, normalized to the total weight in /, J dioI^Lu) is plot- 
ted in Fig. 2 as a function of ^ for different values of 
the disorder strength Vq. We see that the factor of 
in Eq |4^ means that as criticality is approached, almost 
all of the weight in the droplet probability distribution is 
in droplets which do not tunnel. 

For the droplets which are able to tunnel at frequency 
Lj, we find from Eqs. p3 and p6| that 



/l=2(y) 



Co 



c^^y''a{y)C2+ln{J) 



(50) 



Note that in contrast the expression for / in the z = 1- 
case shown in Eq. in the z = 2 case, / does not 
diverge as a; — > 0. As in the z = I case considered above, 
one obtains an expression for /rcst(ti^) by substituting the 
result for / into Eq yielding 



t(^) 



dyy 



-1/2, 



{y^a{y)) 



(51) 



The resulting expression is to good accuracy proportional 
to l/w times logarithms. The physics is that even the 
average of droplets which are able to tunnel is dominated 
by those droplets on the verge of freezing, leading again 
to a superparamagnetic contribution to the susceptibility. 

We have numerically evaluated the integral in Eq. M 
for parameters such that Iq is not too large. Sample 
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results are shown in Fig. 3, which plots the quantity 
Jrest — direst for a relatively small value of the damp- 
ing. The frequency dependence is a consequence of the 
logarithmic factors in Eq the non- vanishing intercept 
as w ^ means that up to logarithmic corrections the 
contribution to the susceptibility arising from this term 
is - 1/T. 

10 I ' , ' , ' , ' , ' 1 




0.02 0.04 0.06 0.08 0.1 

CO/COq 

FIG. 3: Contribution Jrest{(^) ~ ujIrest{oj) (Eq ^l|) of tun- 
nelling droplets to susceptibility integral, plotted vs frequency 
for non-universal constants Vo = 0.5, C2 = 1 and relatively 
small value of damping coefBcient c.^ = 0.1 at ^ = 5 (top 
curve) and ^ = 20 (bottom curve). 



IV. CONCLUSION 

This paper presents an investigation of the possibility 
of quantum Griffiths effects in three dimensional metal- 
lic system which, when 'pure' (non-disordered) is near 
an antiferromagnetic quantum critical point with Ising 
symmetry. For comparison we present also a parallel in- 
vestigation of quantum Griffiths effects in a model of an 
insulating system near a similar critical point. The key 
feature of metallic systems is the dissipative dynamics 
arising from the particle- hole continuum of electrons; in 
the model insulating system the dynamics arc undamped. 
Comparison of the two calculations shows that dissi- 
pation suppresses quantum Griffiths effects completely, 
leaving instead an effectively superparamagnetic behav- 
ior. 

A simple precis of our results follows. Quantum Grif- 
fiths effects are a consequence of randomness: essentially, 
in a random system which is on average in the paramag- 
netic phase, regions ('droplets') may occur in which the 
randomness pushes the system locally to the ordered side 
of the phase diagram, so that local formation of an or- 
der parameter isrfavored. In certain circumstances (first 
noted by McCoya) these droplets may dominate the re- 
sponse. In this situation one may approximately write 



the susceptibility, X; as an average over droplets times a 
susceptibiltiy for each droplet, i.e. 

X = / P(droplet) Xdropict (52) 

J droplets 

We have used simple extremal statistics arguments (sim- 
ilar to those used by Thill and HuseEl) to estimate the 
droplet probability distribution P(droplet) and an ex- 
tension p£ earlier work which studied a particular class of 
dropletsllj to obtain the susceptibility Xdropict of a given 
droplet. We were then able to perform the average over 
droplets and obtain an estimate for the susceptibility. 

This method reproduces the essential features of the 
standard results for quantum Griffiths effects in un- 
damped (insulating) systems, namely that the low T be- 
havior of the susceptibility is governed by a new exponent 
dcs given by the product of the inverse correlation length 
and inverse mean square disorder amplitude V^~^ and 
a non-universal number (which we estimate for the par- 
ticular model we consider). A divergent susceptibility 
results when dcS becomes less than unity, and the results 
are functions only of £^Vq . We note one additional inter- 
esting finding. The standard arguments which produce 
the standard quantum Griffiths results are based on a 
picture of dilute 'droplets' and apply only if the ^ is not 
too large (otherwise the droplets percolate, and an iso- 
lated droplets picture fails). For the model we consider 
we obtain an estimate for the critical value of ^, and find 
that depending on the value of the non-universal factor in 
dcS, droplets may reach the percolation point before the 
Griffiths exponent drops below unity. In other words, in 
the models we consider the existence of a quantum Grif- 
fiths regime (which one may somewhat imprecisely define 
as a divergent susceptibilty arising from fluctuations of 
isolated droplets) is not guaranteed-it may or may not 
occur depending on the value of a non-universal coeffi- 
cient. Sufficiently near a critical point a regime of diver- 
gent susceptibility does of course occur, but the proper 
theory of this regime would have to go beyond the model 
of isolated droplets and treat correctly the mixing of crit- 
ical and griffiths singularities. 

We also found that for undamped systems near antifer- 
romagnetic critical points the amplitude of the divergent 
term in the susceptibility vanishes rapidly as criticality is 
approached, indeed as essentially because each rele- 
vant droplet has a magnetic moment of the order of unity 
and as criticality is approached the droplets get larger in 
size but fewer in number. 

The main new result of our work, however, pertains to 
metallic systems with overdamped (dissipative) dynam- 
ics. For these systems (i.e. for quantum critical phe- 
nomena in metals) the answer is entirely different. The 
physics in the undamped case is a balance between the 
probability of a droplet occurring (which vanishes rapidly 
as the droplet size or amplitude increases) and Xdropict, 
which is of the order of the inverse of the quantum tun- 
nelling rate of the droplet and diverges rapidly as the 
droplet size or amplitude increases. The effect of dis- 
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sipation is to strongly decrease the tunnelling rate, and 
indeed to drive it to zero for droplets larger than a partic- 
ular, amplitude-dependent, size. For relevant parameters 
we find that a non- vanishing density of droplets does not 
tunnel at T = 0; these give rise to a 'superparamag- 
netic' (x ^ susceptibility rather than a quantum 
griffiths (continuously varying exponent) behavior. For 
those droplets which do behave quantum mechanically, 
the effect of dissipation on the tunneling rate is found 
to change the balance between probability and Xdropict 
dramatically. We find that even considering only the 
droplets which can tunnel quantum mechanically, those 
which dominate the integral for x ^-re those which are 
right on the edge of classical (non-tunncUing) behav- 
ior, leading again to superparamagnetism rather than to 
quantum Grifhths behavior. We also find that the depen- 
dence on parameters is different: in the undamped case, 
apart from prefactors the mean square disorder strength 
Vq and the correlation length enter via the combination 
^Vq . In the damped case additional factors of ^ occur 
which drive the system more rapidly to classical behav- 



Our results raise questions about the claimsoEJIi 
that quantum Griffiths effects are important in heavy 
fermion materials, which are precisely three dimensional 
metals with Ising symmetry, typically Jiear antiferromag- 
netic quantum critical points. RefJa contains a phe- 
nomcnological description of data. If the theoretical re- 
sults presented here are accepted, then these data require 
a different, non-Grifhths interpretation. Ref. 113 argued 
that a disordered system near a quantum critical point 
could be mapped onto the disspationless Ising model irua 
transverse field; the results of the present paper and oft3 
indicate on the contrary that dissipation is essential. 

RefJli uses-a novel variant of a technique introduced 
by DotsenkcO to study essentially the same model as 
is studied here. A rather different result was obtained, 
namely that quantum Griffiths effects can be important 
in a reasonable range of the phase diagram even in the 
metallic case. We outline the differences between the 
results found here and_those of Ref.t3. The method in- 
troduced by Dotsenkollj and used by Reitl3 begins from 
a classical theory defined by a functional integral with 
action given by the static term in Eq. |^ and evaluates 
the disorder-swwage by the replica method. Whereas 
other workerglS'llS then used the replicated field theory 
to derive scaling equations for variables including the 
mean disorder strength, Dotsenko argued that one should 
look for spatially localized energetically unstable config- 
urations of the replicated field theory, which correspond 
to local maxima of the replicated action and are to be 
identified with the 'droplets' discussed above. Dotsenko 
shows that the leading nonanalytic contribution to the 
free energy in the vicinity of an assumed T > critical 
point comes from droplets with size of the order of the 
magnetic correlation length, ^; we refer to these hence- 
forth as 'typical droplets'. The authors of Ref.li3 assume 
that the T ^ limit of this classical theory may be 



straightforwardly taken, and then add to this theory es- 
timates of the cl#aamics of 'typical' droplets. The results 
reported in Ref Jl3 disagree in a number of specific details 
with the results presented here, including for example the 
way in which the bare tunnelling rate is estimated. The 
most important differepee, however, is in the interpreta- 
tion of the results. RefJla argues that one should identify 
the boundary of the Grifhths region with the value of 
at which a 'typical droplet' ceases to tunnel. Our analy- 
sis, which involves averaging over all droplets, indicates 
that independent of whether the 'typical droplet' (how- 
ever defined) may tunnel, the susceptibility is dominated 
by droplets which are at or beyond the edge of ceasing 
to tunnel; these give an essentially 'superparamagnetic' 
(x ~ 1/^) behavior, instead of the continuously varying 
exponent characteristic of quantum Griffiths behavior. 

Ref.0 presented a detailed analysis of a different model 
in which spins are added to a pure system which itself is 
far from any critical point. In this model the phase tran- 
sition is disorder-driven: it occurs when the density of 
added spins is high enough that these order; whereas our 
interest here has been in models in which even the non- 
disordered systemia near a critical point. Further in the 
model studied in til the way the disorder is introduced 
means that the local spin amplitude 0o (cf our Eq. |9|) is 
always of order unity, whereas in our treatment the local 
spin amplitude may be considerably smaller. An approx- 
imate mapping between the model considered irU and 
the one considered here may be obtained by setting our 
parameters 0o a-nd ^ equal to unity and considering the 
behavior as the disorder strength Vq is increased (whereas 
we consider a fixed Vq and study the behavior as ^ is in- 
creased). 



r, in a broad qualitative 
Lil are similar to those ob- 
II states that at sufficiently 



Although specific details d 
sense results obtained in Rci 
tained here. In particular, Re: 
low temperature dissipation wiU_suppress the quantum 
griffiths behavior. However, ReftiJ argued that in an ex- 
tremely wide temperature regime could exist in which 
behavior characteristic of the undamped system occurs, 
whereas in the model we consider, for any reasonable 
parameters there|ia no such temperature regime. A cru- 
cial point is thatO focussed on model parameters such 
that the damping coefficient was extremely weak (i.e in 
our notations (see below Eq. H^) they took Cj « 1). 
In this limit, it is plausible that there is a temperature 
regime in which behavior characteristic of the undamped 
model may occur, before finally a crossover occurs to a 
regime (similar to the one we considered) in which damp- 
ing is important. Important avenues for future investiga- 
tion include more detailed studies of the crossovers be- 
tween the weak-damping and order unity damping cases 
aiid between the disorder driven-criticality effects studied 
ir£d and the pure system criticality-driven effects studied 
here, as well as determination of the damping coefficient 
values appropriate to the heavy fermion materials of in- 
terest. 

Our work has the following implications for experi- 
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ment. First, the canonical quantum Griffiths effects are 
due to weak disorder added to a pure critical point. We 
have shown that in the limit of weak disordfitjaiid a pure 
critical point described by the Hertz theorycjO, the dis- 
sipation characteristic of metallic systems changes the 
quantum Griffiths singularities into a kind of superpara- 
magnetic behavior. In other words, as a matter of prin- 
ciple the canonically defined quantum Griffiths behavior 
should not be observable in metals near magnetic quan- 
tum critical points. This suggests that claimsO to have 
observed quantum Griffiths behavior in heavy fermion 
systems should be treated with caution (at least for sys- 
tems with Ising symmetry). Further, we showed that 
the droplets which dominate the susceptibility can tun- 
nel only when the system is not close to criticality, and 
in these cases the droplet size is not much larger than 
the basic scale of the theory. Thus if the susceptibility 
is dominated by the tunnelling of droplets, the picture 
which emerges is more similar to the 'Kondo disorder' 
picture olE3a than it is to the conventionally-defined 
quantum Griffiths picture. Indeed, the experimental 
claims involve 'heavy fermion' systems where the inter- 
action which favors a non-magnfttip phase is the Kondo 
effect. As noted by many authorsE3, the fact that Kondo 
temperatures are exponentially sensitive to system pa- 
rameters means that a slight variation in system param- 
eters can lead to a wide variation in Kondo temperatures. 
The can/jM|Ucal assumption of weak disorder which we and 
otherstlclall3 have made may not be valid for these sys- 



tems. The interplay between quantum criticality and a 
broad distribution of disorder should be treatable by the 
methods introduced here, and seems worth examining. 

A second point is that the very slow dynamics of the 
droplets makes it much easier for them to order. Further, 
in a metallic system the-|droplet-droplet interactions are 
of long range (see, e.g.Efl for a discussion in the context 
of the two dimensional metal insulator transition). For 
this reason we expect that in the presence of disorder the 
actual phase transition at which long ranged order sets 
in is an essentially classical affair, in which droplets lock 
together when the temperature becomes lower than some 
droplet-droplet coupling. 

A third point, perhaps relevant beyond the present 
context, is that (as seen for example in Eq^) dissipation 
can have a crucial effect on bare tunnelling rates: in the 
metallic problem we considered the cruical impediment 
to tunnelling of a droplet was found to be the viscosity 
of the medium, not the energy barrier which had to be 
surmounted. 
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